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1  INTRODUCTION 


There  are  numerous  commercial  codes  available  that  numerically  solve  electric-field  integral 
equations  (EFIE)  for  the  current  and  scattered  fields  produced  by  an  electromagnetic  wave 
incident  on  a  perfectly  conducting  (PEC)  scatterer.  For  research  purposes  and  for  certain 
specialized  geometries,  however,  it  is  often  advantageous  to  have  the  flexibility  and  control 
inherent  in  writing  one’s  own  computer  code.  By  far,  the  simplest  approach  for  numerically 
solving  integral  equations  is  to  use  the  method  of  moments  (MOM)  with  pulse  basis  functions 
and  point  matching.  Unfortunately,  the  EFIE,  unlike  the  magnetic-field  integral  equation 
(MFIE),  cannot  be  solved  accurately  using  pulse  basis  functions  and  point  matching  (as 
demonstrated  in  Figure  5  below),  and  thus  considerably  more  effort  is  required  to  write  a 
computer  code  for  numerically  solving  the  EFIE.  (It  can  be  shown  that  it  is  the  higher-order 
singularity  of  the  EFIE  kernel  that  prevents  an  accurate  solution  using  pulse  basis  functions 
and  point  matching.)  For  open  scatterers  or  many  thin  bodies  the  solution  to  the  MFIE  is 
inderminant  or  unstable,  respectively,  and  it  becomes  necessary  to  use  the  EFIE.  Moreover, 
because  confidence  in  numerical  solutions  is  greatly  enhanced  by  having  two  independent 
numerical  solutions  that  agree  to  within  a  certain  accuracy,  it  is  often  highly  desirable  to 
obtain  the  solution  to  the  EFIE  even  if  the  MFIE  is  also  applicable. 

The  main  purpose  of  this  report  is  to  derive  a  low-order  singularity  electric-field  integral 
equation  (LEFIE)  that  can  be  accurately  solved  using  the  MOM  with  pulse  basis  functions 
and  point  matching.  This  LEFIE,  whose  kernel,  like  that  of  the  MFIE,  has  no  singularity 
greater  than  that  of  the  free-space  Green’s  function,  is  solved  numerically  for  perfectly  con¬ 
ducting  bodies  of  revolution  (BORs)  using  pulse  basis  fucntions  and  point  matching.  Deriv¬ 
atives  of  the  current  are  approximated  with  finite  differences  using  a  quadratic  Lagrangian 
interpolation  polynomial.  This  simple  solution  of  the  LEFIE  is  contingent,  however,  upon 
the  vanishing  of  a  line  integral  that  appears  when  the  original  EFIE  is  transformed  to  obtain 
the  LEFIE.  This  requirement  generally  restricts  the  simple  applicability  of  the  LEFIE  to 
smooth  closed  surfaces.  Bistatic  scattering  calculations  performed  for  scattering  of  a  plane 
wave  by  a  prolate  spheroid  demonstrate  that  numerical  results  comparable  in  accuracy  to 
the  conventional  EFIE  can  be  obtained  with  the  LEFIE  using  pulse  basis  functions  and 
point  matching,  provided  a  higher  density  of  points  is  used  close  to  the  ends  of  the  BOR 
generating  curve  to  compensate  for  the  use  of  one-sided  finite  difference  approximations  of 
the  first  and  second  derivatives  of  the  current. 

The  organization  of  the  report  is  as  follows.  Section  2  contains  the  analysis  of  the  low- 
order  singularity  electric-field  integral  equation  (LEFIE)  solution  of  the  body  of  revolution 
(BOR)  scattering  problem.  It  is  divided  into  several  subsections  beginning  with  the  deriva¬ 
tion  of  the  general  LEFIE  in  Section  2.1  and  its  restatement  for  a  BOR  in  Section  2.2.  The 
solution  of  the  LEFIE  for  a  closed  BOR  using  pulse  basis  functions  and  point  matching  is 
outlined  in  Section  2.3.  Detailed  expressions  for  the  elements  of  the  Z  matrices  that  multiply 
the  column  vectors  of  the  surface  current  expansion  function  coefficients  to  be  determined 
are  derived  in  Section  2.4.  The  Z  matrices  are  treated  by  expressing  them  as  the  sum  of  four 
submatrices  corresponding  to  each  of  the  four  terms  under  the  integral  sign  in  the  LEFIE. 
In  Section  2.5  we  obtain  detailed  expressions  for  the  elements  of  the  V  column  vectors  in 
the  right-hand  side  of  the  matrix  equation  formulation  of  the  LEFIE.  In  Section  2.6  we  ob¬ 
tain  expressions  for  the  currents  induced  on  the  surface  of  a  BOR  by  a  transverse  electric 
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(TE)  and  transverse  magnetic  (TM)  linearly  polarized  plane  wave  in  terms  of  the  solution 
to  the  LEFIE  matrix  equation,  and  in  Section  2.7  expressions  for  the  components  of  the  far 
scattered  field  are  derived.  Section  2.8  discusses  the  choice  of  the  number  of  Fourier  modes 
that  need  to  be  used  in  the  calculations.  Section  2.9  treats  the  approximation  of  the  current 
derivatives  that  appear  in  the  formulation  of  the  LEFIE. 

Section  3  contains  numerical  results  of  calculations  performed  with  a  computer  program 
written  to  implement  and  validate  the  solution  of  the  LEFIE  given  in  Section  2.  A  report 
summary  is  given  in  Section  4. 


2  ANALYSIS 


2.1  Derivation  of  the  Low-Order  Singularity  Electric-Field 
Integral  Equation 

To  derive  the  low-order  singularity  electric-field  integral  equation  (LEFIE)  we  first  derive 
the  ordinary  electric-field  integral  equation  (EFIE)  [1].  On  the  surface  S  of  a  PEC  scatterer 
the  total  tangential  electric  field  vanishes.  The  total  field  is  expressed  as  the  sum  of  the 
incident  field  and  the  scattered  field,  so  that 


n(r)  x  Etot(r)  =  n(r)  x  [Eiflc(r)  +  E8C(r)]  =  0,  r  on  S  (1) 


where  Etot(r)  and  Esc(r)  are  the  total  and  scattered  electric  fields,  respectively,  n  is  the  unit 
normal  vector  to  the  surface  S  at  the  position  r,  and  n  is  assumed  directed  outward  from 
S.  The  scattered  field  can  be  expressed  in  terms  of  a  vector  potential  A(r)  and  a  scalar 


potential  <3>(r)  by 
where 


and 


Esc(r)  =  —juA(r)  —  V4»(r) 

(2) 

A(r)  =  IM)  J  K(r')G(r,  r')dS' 
s 

(3) 

$(r)  =  I  f 

60  { 

(4) 

Here  G(  r,  r')  is  the  free-space  Green’s  function  for  harmonic  time  dependence  exp  (jut)  with 
the  frequency  u>  >  0.  That  is 


G(r,r')  =  G(|r-r'|)  = 


exp(—  jk  |  r  -  r7 1) 
47r  |  r  —  r'  | 


(5) 


where  r  and  r'  are  the  vectors  to  the  field  and  source  points  respectively,  K(r')  is  the  electric 
current  on  S  to  be  determined,  Hq  and  €o  are  the  permeability  and  permittivity  of  free  space 
respectively,  k  —  uj / c  with  c  the  speed  of  light  in  free  space,  and  o  is  the  surface  charge 
density  given  by 

<r(r')  =  •  K(r') .  (6) 
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The  operator  V^-  is  the  surface  divergence  [2,  Appendix  2,  18.].  Combining  (1)  and  (2), 
substituting  (3),  (4)  and  (6),  letting  the  observation  point  r  approach  the  surface  S,  and 
dividing  by  the  free-space  impedance  Z0  =  (A*o/e o)1^2  we  obtain 

^n(r)x^[*:2K(r')G(r,r')-(Vi-K(r'))V'G(r,r')]<iS'=^xEi’“(r).  (7) 

s 

The  “o”  on  the  integral  sign  in  (7)  indicates  that  a  small  “principal  area”  isolates  the 
singularity  of  the  Green’s  function  from  the  surface  integration.  The  form  of  the  EFIE  in 
(7)  is  conditional  upon  the  choice  of  the  principal  area  being  a  circle  with  the  singular  point 
at  its  center  (or  another  principal  area  that  is  adequately  symmetric  with  respect  to  the 
singular  point)  [3].  Equation  (7)  is  the  conventional  form  of  the  EFIE  for  the  current  K  on 
the  surface  S  of  a  PEC  scatterer. 

The  gradient  operator  acting  on  the  free-space  Green’s  function  in  (7)  results  in  a  higher 
order  singularity,  1/|  r  —  r'  |2,  instead  of  the  1/|  r  —  r'  |  singularity  for  the  free-space  Green’s 
function  itself.  What  we  want  to  do  is  to  recast  the  conventional  form  of  the  EFIE  into  a 
new  form  that  has  no  singularity  higher  than  that  of  the  free-space  Green’s  function.  To  do 
this  we  begin  by  writing  the  gradient  operating  on  the  free-space  Green’s  function  as  the 
sum  of  the  surface  gradient  V5  and  the  gradient  in  the  normal  direction  so  that  (7)  becomes 


|n(r)  x  f^K(r')G(r,  r')  -  (V's  •  K(r'))V'sG(r,  r')  - 


s 


( V's  ■  ^  X  E‘~(r).  (8) 

The  term  in  (8)  with  the  normal  derivative  of  G  includes  the  cross  product  n(r)  x  n(r')  and 
does  not  have  a  singularity  higher  than  G  itself  at  r  =  r'.  Focusing  on  the  surface  gradient 
term,  we  use  the  vector  identity  [2,  Appendix  2,  26.] 


Vs(AB)  =  AVSB  +  BVSA 


(9) 


to  write 

(Vs  ■  K )V'sG  =  V'S[(V^  •  K)G]  -  GV's(V's  ■  K). 
We  then  use  a  Gauss  integral  theorem  [2,  Appendix  2,  43.]  to  obtain 


(10) 


/  Vg[(Vg  •  K)G]dS'  =  -j  J'(V'S  •  K)Gn'dS'  +  J (Vs  •  K)Gm,|dc'|.  (11) 


In  (11) 


J' 


R[  + 


(12) 


where  R[  and  R^  are  the  principal  radii  of  curvature  at  the  surface  point  r',  and  C  is  a  set 
of  closed  craves  (defined  below)  on  the  surface  of  the  PEC  scatterer.  The  unit  vector  m', 
defined  at  each  point  on  the  curves  comprising  C,  is  in  the  tangent  plane  to  the  surface  at 
the  point  and  perpendicular  to  the  curve.  Primes  are  used  to  denote  dependence  on  the 
source  point  r'. 


3 


Since  (11)  is  the  key  step  in  obtaining  a  new  form  of  the  EFIE  it  is  important  to  under¬ 
stand  the  meaning  of  S  and  C  in  applying  it.  We  assume  that  the  surface  S'  of  a  general 
PEC  scatterer  can  be  divided  into  a  finite  number  of  open  subsurfaces  sharing  bounding 
curves  in  common  with  one  another,  such  that  on  each  subsurface  the  surface  charge  density 
<7(r')  =  -l/(ju)Vs  •  K(r'),  its  surface  derivative,  and  J'  are  continuous  and  integrable,  and 
n'  and  m'  are  continuous.  In  applying  (11)  to  the  surface  S  of  the  PEC  scatterer,  S  is  to  be 
regarded  as  the  superposition  of  these  subsurfaces,  and  the  surface  and  line  integrations  per¬ 
formed  separately  for  each  subsurface  with  its  bounding  curve.  The  set  of  bounding  curves 
comprise  C.  The  unit  vector  m'  on  each  of  the  bounding  curves  comprising  C  points  away 
from  the  subsurface  it  encloses.  The  surface  of  a  closed  finite  cylinder,  for  example,  is  to  be 
regarded  as  the  superposition  of  three  subsurfaces,  the  side  cylindrical  surface  and  the  two 
disks  at  either  end.  The  surface  of  a  sphere,  a  simple  smooth  scatterer,  can  be  regarded  as 
the  superposition  of  two  hemispheres. 

Combining  (8)  and  (11)  we  obtain  a  new  form  of  the  EFIE  with  a  self-term  singularity 
equal  to  that  of  the  ffee-space  Green’s  function 

|n(r)  x  /{[lfc2K(r')  +  VS(V'S  •  K(r'))  +  /(V's  •  K(r'))n(r')]G(r,  r1) 


(V's  •  K(r'))a6^’,r-n(r')}<i,S'  -  |n(r)  x  J (V's  ■  K(r'))G(r,r')m(r')|<fc'| 


n(r)  x  E”(r)  (13) 

Zo 

which,  unlike  the  original  EFIE  in  (7),  is  not  conditional  upon  the  shape  of  the  principal 
area  used  to  isolate  the  singularity  of  the  Green’s  function. 

Now  our  recasting  of  the  conventional  EFIE  (7)  into  this  new  low-order  singularity  form 
(13)  is  motivated  by  the  desire  to  be  able  to  solve  the  EFIE  with  the  method  of  moments 
using  pulse  basis  functions  and  point  matching.  Therefore  we  would  like  to  avoid  the  line 
integral  over  C  in  (13).  Assume  that  the  surface  S  of  the  scatterer  is  closed  and  smooth 
enough  so  that  a(r')  =  -l/{ju)Vs  •  K(r'),  its  surface  derivative,  and  J'  are  continuous 
and  bounded,  and  n'  and  m'  are  continuous  over  S.  Then  (13)  can  be  applied  to  any  two 
contiguous  open  surfaces  (with  a  smooth  bounding  curve)  comprising  the  closed  surface,  the 
two  line  integrals  over  the  common  bounding  curve  of  the  two  contiguous  surfaces  cancel, 
and  (13)  reduces  to 


x  j{[k2 K(r')  +  V'S(V^  •  K(r'))  +  J'(VS  •  K(r'))n(r')]G(r,  r') 

5 


— (V5  •  K(r')) 


dG{  r,r') 
dri 


n(r  ')}dS' 


n(r)  x  Einc(r) 
£0 


(14) 


If  the  scatterer  has  an  open  surface,  edges,  tips,  or  any  boundary  where  a  or  J'  become 
singular  or  discontinuous,  or  n'  and  m'  become  discontinuous,  the  line  integral  over  C  cannot, 
in  general,  be  omitted  and  the  numerical  solution  to  the  low-order  singularity  electric-field 
integral  equation  (LEFIE)  (13)  may  be  more  complicated  than  that  of  the  original  EFIE. 
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Accordingly  we  will  restrict  our  treatment  of  the  LEFIE  in  this  report  to  PEC  scatterers 
with  closed  smooth  surfaces  such  as  spheroids  for  which  <j,  its  surface  derivative,  and  J'  are 
continuous  and  bounded,  and  n'  and  m'  are  continuous,  so  that  the  line  integral  over  C  can 
be  omitted  and  the  LEFIE  takes  the  form  (14)- 


2.2  LEFIE  for  a  BOR 

We  seek  to  determine  the  surface  current  and  the  far  scattered  field  of  a  perfectly  electrically 
conducting  (PEC)  closed  body  of  revolution  (BOR)  excited  by  an  incident  plane  wave. 
The  geometry  of  the  BOR  is  shown  in  Figure  1.  Circular  cylinder  coordinates  (p,  0,  z )  are 
employed  with  (p,  0,  z)  denoting  the  corresponding  unit  vectors,  and  with  the  z  axis  chosen 
as  the  axis  of  revolution.  The  origin  of  the  circular  cylindrical  coordinate  system  lies  on 
the  z  axis  but  does  not  necessarily  coincide  with  the  lower  pole  of  the  BOR  as  in  Figure 
1.  The  coordinates  (i,0),  with  t  the  path  length  along  the  generating  curve  of  the  BOR 
from  the  lower  pole,  form  an  orthogonal  curvilinear  system  on  the  surface  S  of  the  BOR;  the 
corresponding  unit  vectors  are  (t,  0) .  Figure  2  shows  the  propagation  vector  kmc  =  fcktnc  of 
the  incident  plane  wave.  The  propagation  vector  is  assumed  to  lie  in  the  xz  plane  (0  =  0), 
with  — kmc  making  an  angle  of  9tnc  with  the  positive  z  axis  and  with  kxinc  <  0  so  that 

kinc  =  — fc(sin  0incx  +  cos  9incz) .  (15) 

Also  shown  in  Figure  2  are  the  spherical  polar  angles  of  the  far  field  observation  point 
r/ar  _  ^  p/ar)  0/or)  and  the  associated  unit  vectors  0^  and  iff  .  For  TM  illumination  the 
incident  electric  field  is  given  by 

Einc  =  kZ0exp(-jkinc  •  r)0m°  (16) 


while  for  TE  illumination 

Einc  =  kZ0exp(-jkinc  •  r)#™.  (17) 

In  (16)  and  (17)  r  is  the  vector  from  the  origin  to  any  point  in  space  and  the  factor  of  kZ0 
is  inserted  to  simplify  the  analysis. 

For  a  BOR  (14)  can  be  replaced  by  the  equivalent  pair  of  equations 


t(r)  •  /{[fc2K(r')  +  V'S(V'S  ■  K(r'))  +  J'(V'S  ■  K(r'))n(r')]G(r,r') 


and 


|«(r)  •  j{\e. K(r')  +  V's(Vi  •  K(r-))  +  J\VS  ■  K(r'))fl(r')]G(r,r') 


(18a) 


-(V..K  M)^(rW  =  M_p> 


(18b) 
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The  LEFIE  for  a  PEC  scatterer  with  a  closed  smooth  surface  in  (14)  and  for  closed 
smooth  BORs  in  (18),  like  the  original  EFIE  in  (7),  produce  a  unique  solution  for  the  surface 
current  K  except  at  frequencies  equal  to  the  resonant  frequencies  of  the  cavity  formed  by 
the  closed  surface  S  of  the  scatterer.  These  spurious  resonances  can  be  eliminated  from  the 
LEFIE  in  the  same  way  they  have  been  eliminated  from  the  EFIE.  For  example,  the  LEFIE 
can  be  combined  with  the  magnetic-field  integral  equation  [4]  or  added  to  a  corresponding 
LEFIE  that  is  satisfied  on  a  dual  surface  just  inside  the  surface  S  of  the  scatterer  [5]-[9].  To 
concentrate  on  the  subject  of  lowering  the  order  of  the  singularity  of  the  EFIE  and  not  on  the 
details  of  these  methods  for  eliminating  spurious  resonances,  we  shall  choose  frequencies  in 
our  numerical  examples  that  are  sufficiently  far  from  any  cavity  resonance  to  avoid  numerical 
instabilities. 

2.3  Solution  of  the  LEFIE  by  the  Method  of  Moments 

To  solve  the  LEFIE  (18)  for  the  surface  current  K  we  begin  by  expanding  K  in  a  Fourier 
series 

K(fV)=  £  [K'Smt'A')  +  K*(?)i>(#)}  (19) 

n=—N 

The  choice  of  the  value  of  N  is  discussed  in  Section  2.8.  To  obtain  separate  integral  equations 
for  each  of  the  Fourier  modes  we  multiply  both  sides  of  (18a)  and  (18b)  by  e  Jm<^,ra  = 
0,  ±1,  ±2,  •  •  • ,  and  integrate  with  respect  to  4>  from  —tv  to  7r.  As  will  be  seen  below,  the 
integrands  of  the  LHSs  of  (18a)  and  (18b)  are  of  the  form  F(t,t  ,<j>  —  <f>)e?  ^ ,  if^the  dot 
products  of  t  and  <f>  are  taken  inside  the  integral  signs.  Noting  that  dS'  =  p'dt'dfi  and 
performing  the  integration  with  respect  to  <ft'  from  — 7T  to  7r  as  well  as  the  integration  with 
respect  to  <f>  we  then  have 

j  d(f)e~jm<l>  j  d0V'n^'F(t,  t',  <£'-<£)  =  j  d(pej(n-rn)4'  j  dcjJeP^'^F^,  t',  <}>'  -  4>) 

— 7T  —  7T  ~7r  ~~n 

7T 

-  2t rSnm  J  d<f>'^'F(t,  t',  <t>')  (20) 

—7 r 

where  the  Kronecker  delta  6nm  equals  0  for  m  n  and  equals  1  for  m  —  u.  The  orthogonality 
of  the  Fourier  modes  thus  enables  separate  integral  equations  to  be  obtained  for  each  Fourier 
mode  ein<t>. 

Following  the  (j)  and  (j)'  integrations  we  approximate  K^(t')  and  Kn{t )  by  pulse  basis 
functions  defined  as  follows.  The  generating  curve  of  the  BOR  is  parametrized  in  terms  of 
t,  the  distance  along  the  curve  measured  from  the  lower  pole  of  the  BOR.  For  each  value  of 
t,  the  corresponding  point  on  the  generating  curve  is  given  by  [p(t),z(t)\.  A  set  of  M  +  1 
points  p\,P2,  •  •  •  >Pm+i>  is  chosen  to  discretize  the  generating  curve  with  p*  =  (p\,Zi)  the 
lower  pole  of  the  generating  curve  corresponding  to  t  =  0,  p*M+i  =  (Pm+d  zm+i)  the  upper 
pole  of  the  generating  curve,  and  with  t*  >  t*_v  The  generating  curve  is  approximated  by 
straight  line  segments  between  adjacent  points.  The  midpoints  of  the  approximating  straight 
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line  segments  axe  given  by 


(Pi.Z.)  =  ( 


Pl  +  fjg  4+ik 


)•*“ 


1,2, •••,*/ 


(21) 


and  the  length  of  the  ith  straight  line  segment  denoted  by 


di  = 


(Pi+i  _  Pi)  +  (zi+ 1  -  zt) 


!*1  1/2 


(22) 


For  calculation  purposes  the  discretized  generating  curve  then  completely  replaces  the  origi¬ 
nal  generating  curve  and  the  parameter  t  now  becomes  the  length  along  the  discretized  curve 
from  the  lower  pole  instead  of  the  length  along  the  original  generating  curve  from  the  lower 
pole.  Thus,  for  example,  t(pl)  =  0,  t{p%)  =  dx,  t(pz)  -  dy  +  d2,  etc.  A  pulse  basis  function 
Pi(t)  is  defined  as 

Pi(t)  = 

Then 

M 

KWkZKVMV)  (24a) 


( 


o,  t<n,t>n+ 1 
i,  n<t<t*+1  . 


(23) 


and 


M 


K*(t)  *Y,K*M)Pin 
*= 1 


(24b) 


A  set  of  2 M  equations  for  the  2 M  unknowns  Kln(ti),  K%(ti),  i  =  1, 2,  •  •  • ,  M,  is  then  obtained 
by  using  point-matching:  the  LHS  and  RHS  of  the  integral  equations  for  each  Fourier  mode 
are  equated  at  t  =  U,  i  =  1, 2,  •  •  • ,  M.  This  set  of  2 M  equations  can  be  expressed  in  matrix 
form  as 

T  [rrtti  {rrtsi  1  r  t rt  1  I  \rt  1 

n  =  0,  ±1,  ±2,  •  •  • ,  ±N.  (25) 


[Z»\  1 

K' 

V*  1 

T  n 

m  J 

Kt  J 

\T<t> 
n  - 

In  (25)  the  [Z^9],  p,q  =  t  or  <j),  are  M  x  M  matrices  obtained  from  the  LHS’s  of  (18a)and 
(18b).  The  index  p  corresponds  to  the  external  dot  product  factor  t  or  </>,  and  the  index  q 
corresponds  to  the  t  or  <f>  component  of  K.  The  ith  row  of  [Z^9]  corresponds  to  the  value  of 
the  observation  point  ti,  and  the  jth  column  of  [Z?q]  corresponds  to  Kq(tj).  The  ith  value 
of  the  M  x  1  vectors  K‘  and  equals  K^U)  and  K£(U)  respectively.  The  vectors  V* 
and  on  the  RHS  of  (25)  contain  the  values  of  the  RHS’s  of  (18a)  and  (18b),  respectively, 
evaluated  at  t  =  ti,  %  =  1, 2,  •  •  • ,  M,  following  multiplication  by  e~9n<t>  and  integration  with 
respect  to  </>  from  —  •n  to  7r. 


2.4  Expressions  for  the  Elements  of  the  Z  Matrices 

Detailed  expressions  will  now  be  obtained  for  the  elements  of  the  Z  matrices.  The  contri¬ 
bution  of  the  dot  product  of  t  or  with  each  of  the  four  terms  under  the  integral  sign  in 
(18a)  and  (18a)  will  be  considered  in  turn.  Some  useful  relationships  will  first  be  established. 
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Letting  v  and  v'  be  the  angles  measured  positive  clockwise  from  the  positive  z  axis  tc 
t*  respectively  (see  Fig.  3),  it  follows  that 

t  =  cos  v(t) z  +  sin  v(t)p  =  sin  v(t)  cos  </>x  +  sin  v(t)  sin  <py  +  cos  v(t)z, 

(j>  =  —  sin  (fyk  +  cos  4>y 

A  f  A  / 

and  similarly  for  t  and  <f>  ;  and  hence  that 


so  that 


cos  v'  cos  0'x  +  cos  v'  sin  4>  y  -  sin  v'z 

t  •  =  sinu(t)  sin  v'{t')  cos{(f>  —  4>)  +  cos  v(t)  cos  v'{t') 

t  •  (j>  =  —  sinu(t)  —  <f>) 

<}>  •  t'  =  sin  v\t')  sin(0'  -  <f>) 

<j>  •  <j>  =  cos(4>'  —  4> ) 

t  •  n'  =  sin  v  cos  v'  cos {<j>'  -  <f>)  -  cos  v  sin  v' 

<(>  •  n'  =  cos  v'  sin (</>'  —  <f>). 

It  is  also  simple  to  show  that 

|  r  -  r'  |  =  [p2  +  pn  -  2 pp'  cos(<j)'  -</>)  +  (z  -  z /)2]  / 


and 


=  [(p  -  p')2  +  (z  ~  z'Y  +  App'  sin2 '  ■ 


Now  let 


IT  =ij  dfrT*”*] p(r)  •  /  fc2^  K{t')e?n+\{r')G{\  r-r'  \)dS' 
k  J  J  n 

— 7T  O 

=  jk  j  d4>e~im* p(r)  •  J p'dt'  j  d<t>'  J2 r-r'\),p,q  =  t  or  <j>. 

— 7T  — 7T  U 

where  the  integration  with  respect  to  t'  is  over  the  length  of  the  generating  curve 
BOR.  Then,  using  (5),(28),(27a)-(27d),  and  (20)  we  obtain 

If  =  jk 2  J dt'p'Ki (*') [sin v sin v'G2,n (p,  p',z  -  z') +  cosv cos v'Ghn (p,  p’,z-z')] 

if  =  - k 2  J  dt'p'Klit')  sinu/G3)„(p,  p',  z  -  z') 
if  =  A:2 sin  v  j dt' p' K*(t')Gz>n{p,  p',  z  -  z') 
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t  and 

(26a) 

(26b) 

(26c) 

(27a) 

(27b) 

(27c) 

(27d) 

(27e) 

(27f) 

(28) 

(29) 
of  the 

(30a) 

(30b) 

(30c) 


z 


Figure  3:  Defining  geometry  for  the  angle  v  in  the  pz  plane. 
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and 


where 


and 


with 


It*  =  jk2 j  dt'p’K£(t')G2tn(p,  p',  z  -  z’) 

t r 

Gitn(p ,  p',  z-z')  =  J  G0(R)  cos {n4>')d(j)' 

0 

7T 

G2>n{p,  p',  z-z')  =  J  Gq{R )  cos (ruff)  cos 4>'d(f)' 

0 

7T 

G3iJp,fi,  z  -  z')  =  J  Go(R)  sm(rup')  sin  fid# 


Go(fl)  = 


exp(-jkR) 

fc# 


(30d) 

(31a) 

(31b) 

(31c) 

(31d) 


and  /?  =  |  r  —  r#  |  given  by  (28).  In  the  arguments  of  the  Gfc,n>  fc  =  1>  2, 3,  the  distances  p 
and  z  are  functions  of  the  observation  point  coordinate  t,  and  p'  and  z'  are  functions  of  the 
source  point  coordinate  t'. 

The  contributions  of  the  I[q  to  the  Z?q  matrices  are  then  (see  Section  2.3) 

{ZZUij  =  j k2dj [sin  Vi  sin  VjG2,n (pi ,  Pj ,  Zi  -  Zj)  +  cos  Vi  cos  VjG1>n(pi,  Pj,  Zi  —  Zj)\  (32a) 

[Znh ;ij  =  ~k2djPj  sin  VjG^n(pu  Pj ,  -  Zj)  (32b) 

lztnh w  =  k2  sin  VidjPjGztn{pu  Pj ,  Zi  -  Zj)  (32c) 

and 

[Zf  ]  1*  =  Pi ,  Zi  -  Zi).  (32d) 

In  (32)  the  subscripts  i  and  j  indicate  observation  point  and  source  point,  respectively,  and 
the  discretized  values  of  sin  u  and  cos  v  are  given  by 


and 


sm  Vi  = 


COS  Vi  = 


P*+i  -  P* 


di 


4+i  ~  4 


(33a) 

(33b) 


Next  let 

I?  =  {}  d<t>e-im* P(r)  -fVsWa-'E 


G(|r-r'|)dS' 


=  J-j  d<t>e~jm<t> p(r)  •  J p'dt'  J  d<t>'Vs  Vs  •  E  Kqn(t'yn+  q(r')  G(|  r  -  r'  |), 


p,  q  —  t  or  <f). 


(34) 
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Making  use  of  the  relationships  [2,  Appendix  2,  73.  and  74.] 


r?'. 


(35a) 

]  = 

(35b) 

(35c) 

& 

1 _ 1 

i'  -  ^K*(t ')£')  e** 

(35d) 

V's  [V8  •  i'Ktifle**]  = 

and  (5),(28),(27a)-(27d),  and  (20),  we  obtain 

I 2  =j  J dt'p'-Jt 7  [sin v  sin v'G2,n(p,  p',  z-z')+  cos v cos v'Gi<n(p,  p',z-  z')] 

./ 1  ^  r  J  r/t  /V 


+jnsinw  J dt'-—\p'Ktn(t,)]G^n(p,  p',  z  -  z') 

if  =  -/ -  *0 
-n  /  tf—KWfe.hp',*-  *0 


(36a) 


(36b) 


=  ~n  J  [Sin  U  S^n  V>^2,n{p,  p',  Z  —  z')  +  COS  U  COS  V,Gfl,n(p,  /t/>  2  “  *01 


-n2  sin  t;  J  dt'Kn^  -G3tn(p,  p',  z  -  z') 


(36c) 


and 


r<£4> 


-jn  J  dt'p‘ 


i  j  d 

dt' 


*m 


sin*/G3in(p,p',z  -  z')  -  jn2  J  dt' ^K*(t')G2,n{p,  p' ,  z  -  z') 


(36d) 

Before  obtaining  the  contribution  of  if9  t°  [Z£9]  ^  ^  necessary  to  consider  the  derivatives 
appearing  in  (36).  The  details  of  the  derivative  approximations  will  be  treated  in  Section 
2.9.  For  our  purposes  here  it  suffices  to  note  that  a  quadratic  Lagrangian  interpolation 
polynomial  is  used  to  obtain  approximations  for  both  the  first  and  second  derivatives.  For 
the  first  derivatives 

d  r  ,rr, ,.,.i 


and 


d_ 

dt' 


1 

X/  Cj+kPj+kKnj+k 

(37a) 

k=-l 

(37b) 

13 


where  we  have  denoted  K^tj)  and  K&tj)  by  Kln>j  and  K*j  respectively.  For  the  second 
derivative 

*  [p'KW)]t.,t  “  E  t38) 


dt' 2 


fe=-i 


Since 


I  ^  W  -  JOT  w  <39) 

£  -  r”-i^  (40) 

where  we  have  used  (33)  to  approximate  dp' /dtf .  In  (37), (38),  and  (40),  when  j  =  1,  k  is 
summed  from  0  to  2,  and  when  j  =  M,  k  is  summed  from  -2  to  0.  Using  the  derivative 
approximations  (37),  and  (40)  we  then  obtain  the  following  contributions  of  the  If  to  the 
Zf  matrices: 

[Znh;i,j+k  =  jdjPj  [sin Vi  sin  VjG2,n(pi, Pj,  zi  ~  zj)  +  cosvi  cos VjGitn(pi, Pj,  z{  -  Zj)] 

•  ( —c"+kpj+k  -  \ sin vjc'j+kpj+kS)  +  jn sin VidjG3,n(pi,  Pj, zi  ~  zj)-kcj+kPj+k  (41a) 
\Pj  Pj  /  ^3 

\Znh;ij+k  =  -djpj  sin  VjGztn{fiii  Pj,  zi  -  Zj)  (jf.cj+kPj+k  ~  ^  smVjc'j+kpj+^j 

TldjG2,n{,Pi,  Pj,  zi  zj)^j+kPj+k  (41b) 

f  1 

[znh-,i,j+k  =  -ndjpj  [sin  Vi  sin  v jG2, n(pi,  Pj,  zi  ~  zj)  +  cos  v{  cos  VjGi,n  (p, ,  pj ,  2*  -  zj)}  cj+fc^“ 


-n2  sin Vidj — G3n (pi ,  pj,  z{  -  Zj)  |fc=0 
Pj 


(41c) 


and 


[zt%;i,j  =  -Jndjpj  sinVjG^n(Pi,  Pj,  zi  -  Zj)  -  jn2-^G2,n(Pi,  Pj,  zi  -  Zj).  (41d) 

Pj 

The  index  k  in  these  expressions  takes  the  values  —1, 0, 1  for  j  =  2, 3,  •  •  • ,  M  —  1,  and  the 
values  0, 1, 2  and  -2,  -1, 0  for  j  =  1  and  j  =  M,  respectively. 

Next  let 

If  =  3-  j  d<pe~jm*p(r)  ■  f  J'VS  ■  q(r')J2Kn(t'yn,i>'  n(r')<?(l  r  “  r'  l)rf5' 

ro  J  J  L  n 


|  j  dfe-** p(r)  ■  J  p'J'tfW  j  d>!>'Vs  ■  q(r')  £ 

_ <rr  —  7T  ' 


n(rOC7(|  r  -  r'  |), 
p,q  =  toT(f).  (42) 
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Using  (35a) ,  (35b) ,  (5) ,  (28) ,  (27e) ,  (27f) ,  and  (20)  we  obtain 

1$  =  j  f  dt' ^-\p' Kl(t')]J{iI)  [sin v cos v'G2,n(p,  p\  z-z')-  cost- sin v'Gltn(p,  p',  z  -  z')] 

d  (43a) 

if  =  -f  df^[p'KtMJ(t')cOSv'G^(p,p',z-z')  (43b) 

I**  —  —n  J dt'K*(t')J(t')  [sin v  cos v'G2,n(p,  p',  z  -  z')  -  cosvsmv'G^n(p,p',  z  -  z')]  (43c) 
and 

1$*  =  ~jn  J  dt'K ’+  (t')J(t')  cos  v'G3<n(p,  p',  z-z').  (43d) 

Using  the  first  derivative  approximation  (37)  the  corresponding  contributions  to  the  Z%q 
matrices  are  then 


[Zn  ]3;i,j+fe  —  jdjj(tj )  [sin  Vi  COS  VjG 2, n{pi,  Pj,  Zj) 

-  cos  Vi  sin  VjGlin(pi,  pj,  z{  -  Zj)]c'j+kpj+k  (44a) 

\^n  ]3;t,j+fc  =  djJ(tj)  COS  VjGstn(pi,  Pj,  Z{  Zj)Cj+kPj- |-fc  (44b) 

[Znh-ij  =  —ndjJ(tj)[sm  cos  vjG2,n(Pi,  Pj,  z{  ~  Zj)  -  cosv,  sin  VjGlin(pi,  pj,  Zi  -  Zj)]  (44c) 

and 

[Z**h -,ij  =  -jndjJ(tj)  cosvjG3<n{pi,  Pj,  Zi  -  Zj).  (44d) 

Finally  let 


.  7T  r 

I?  =  /  #e-«p(r)  •  j  Vs  ■  q(r') 


3G(|r-r'|) 

dri  K  ' 


=  -3-j  p(r)  ■  j  flit  j  d4'V's  ■  [q(r')  £  ~  'W'). 

—7 r  — 7 r 

p,  q  =  t  or  (f>. 

Using  (33),(5),(28),(27e),(27f),(20),  (26c),  and  noting  that 

dG( |  r  -  r'  |) 


dn' 


=  VG( |  r  -  r'  |)  •  n(r') 


V'G(|  r  -  r'  |)  =  i_M+ie-J**(r  —  r'),  R  =  t-  r' 


(45) 

(46) 

(47) 


we  obtain 


/d 

dt'— [p'K^(t')] [psin  v  cos2  v'H^n{p,  p' ,  z-z') -pcos  v  cos  v'  sin  v'H2,n{p,  p' ,  z-z') 

-  sin  vp' cos2  v'H2,n(p,  p',  z-  z')  +  cos  vp'  cos v' sin  v'Hltn(p,  p',  z  -  z') 

— (z  —  z') sin v cost/ sin v' H2in(p, p' , z  —  z')  +  cos v(z  —  z') sin2  v'H\^n(p, p' ,  z  —  z')\  (48a) 
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if  =  k2 1 dt'—ip'K^t'^lpcos2  v'H5>n(p,  p',  z  -  z')  -pcos2  v'H3>n{p,  p',  z  -  z') 

— (z  —  z')  cos  v'  sin  v'H3>n(p,  p'  ,z  —  z')\  (48b) 

if  =  nk 2  J  dt'Kft')  [p sin  v cos2  v'H4<n(p,  p',  z-z')-p cos  v  cos  v'  sin  v'H2>n(p,  p' ,  z  -  z') 

-  sin  vp'  cos2  v'H2,n(p ,  p',  z  -  z')  +  cos  up'  cos  v'  sin  v’Hhn(p,  p' ,  z  -  z') 

— (z  —  z') sin v cos v' sin v'H2<n (p,  p\  z  —  z')  +  cost >(z  —  z )  sin2 v  >n(p,  P  >  ^  —  z  )]  (48c) 

and 

if  =  jk2n  [ dt'Kft') [p cos2  v'H5,n{p,  p',  z  -  *')  -  pcos2  t/tf3,n(p,p',  z  -  z') 

(48d) 
(49a) 
(49b) 
(49c) 
(49d) 
(49e) 

(49f) 

and  R  =  |  r  -  r'  |  given  by  (28).  Using  the  first  derivative  approximation  (37)  the  corre¬ 
sponding  contributions  to  the  Z^q  matrices  are  then 

[Z^U-ij+k  =  -jk2dj[pi  sin  cos2  VjH4in(pi ,  pj,  zt  -  Zj)—pi  cos  u*  cos  Vj  smVjH2>n(pi,  Pj ,  z*  -  Zj) 

—  sin  ViPj  cos2  VjH2<n(pi,  pj,  Zi  —  Zj)  +  cos  ViPj  cos  Vj  sin  VjH^n{pi,  pj,  Zi  —  Zj) 

-( Zi  -  Zj)  sinUi  cos  Vj  sin  VjH2>n(pi,  pj,z{-  Zj) 

-1- cos  Vi{zi  -  Zj)  sin2  VjHUn(pi,  Pj,  z,  -  Zj)]c'j+kPj+k  (50a) 

[ZfU-,i,j+k  =  k2dj[pi  cos2  VjH^n(pi,  pj,  Zi  -  Zj)  -  pi  cos2  VjH3<n(pi,  pj,Zi  -  Zj) 

— (z  -  z')  cos  Vj  sin  VjH^n(pi,  Pj,  Zi  -  Zj)]Cj+kpj+k  (50b) 


In  (48) 


with 


J 

-(z  -  z')  cos v' sin v'H^n(p,  p' ,  z  -  z')\. 

7T 

z')  =  J  H„(R)  cosWW 
0 
7T 

H2i„(p,  p',  z  -  z')  =  J  H0(R)  cos {n<t>')  cos  <fd<f 
0 

7T 

H3,n(p,  p',  Z-  z')  =  J  Ho{R)  sin (n<f>')  sin  <fd<f 
0 
7T 

H4>n(p,  p',  Z~z')  =  J  Ho(R)  cos(n<j/)  cos2  4>'d<p' 

0 

7T 

H5tn(p,  p',  z-z')  =  J  H0(R)  sin (ruf>')  cos  </>'  sin  tfdtf 


Ho (R)  =  eM-jkR) 
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[Zn\4\i,j  =  nk2dj[piSinViCos2  VjH4jn(pi,  pj,  Zi  -  zj)  -  pi  cos  cos  vj  sin  VjH2tn(Pi,Pj,  -  Zj) 
-  sin  ViPj  cos2  VjH2,n(Pi,  Pj,  Zi  —  Zj)  +  cos  ViPj  cos  Vj  sin  VjH1<n{pi,  pj,  Zi  -  Zj) 

-( Zi-Zj )  sin  Vi  cos  Vj  sin  VjH2,n(pi,  Pj,  Zi  -  ^)+cos  Vi(zi  -  Zj)  sin2  VjHhn(pi,  pj,  Zi  -  Zj)]  (50c) 
and 

[Z£%,j  =  jk2ndj[pi  cos2  VjHb>n{pu  pj,  z{  -  zj)  -  pi  cos2  VjH3in{pi,  pj,  z{  -  Zj) 

-( Zi  -  Zj)  cos  Vj  sin  VjH3>n(pi,  pj,  z{  -  Zj).  (50d) 

Equations  (32),  (41),  (44),  and  (50)  give  the  contributions  of  the  four  terms  under  the 
integral  sign  of  (18)  to  the  Z£q  matrices  in  (25)  so  that 


[zr i  =  [zrii  +  [zri2 + wh + wu,  p,? = *  <*  <*>. 

(51) 

It  is  easy  to  see  from  (31),  (49),  (32),  (41),  (44)  and  (50)  that 

lZ—nl  =  ffl] 

(52a) 

[Z?J  =  -[Z*] 

(52b) 

1 

II 

1 1 

(52c) 

and 

ic 

ii 

tf 

(52d) 

2.5  Expressions  for  the  Elements  of  the  V  Vectors 

Now  that  we  have  obtained  detailed  expressions  for  the  elements  of  the  Z  matrices  we  turn 
to  the  elements  of  the  column  vectors  and  in  (25).  Recall  that  the  V  vectors  represent 
the  RHS’s  of  (18)  multiplied  by  exp (—jncp)  and  integrated  with  respect  to  4>  from  —7 r  to  7r 
and  evaluated  at  the  M  points  t  =  U,  i  —  1, 2,  •  •  • ,  M,  and  that  the  incident  electric  field  is 
given  by  (16)  and  (17)  for  TM  and  TE  polarization,  respectively.  Then 


V^  =  ^r  d<f>e-^p(ti,  <f>)  ■  kZ0  exp[kinc  •  r (U,  </>)]qinc,  p  =  t  or  <f>,  q  =  9  or  0.  (53) 

Zo  J 

— 7T 

Using  (26a)  and  (26b)  it  is  simple  to  establish  the  unit  vector  relations 


t  •  0mc  =  cos  9mc  sin  v(t)  cos  (0  —  0inc)  —  sin  6inc  cos  v{t) 

(54a) 

0  •  0inc  =  -  cos  9inc  sin (0  -  0ir*c) 

(54b) 

t  •  0mc  =  sin  v(t)  sin (0  —  0tnc) 

(54c) 

and 

0  •  0inC  =  cos  {(f)  -  <f)inc). 

(54d) 

Also,  using  (15) 

kmc  •  r(ti,  4>)  =  — fc(sin  6tncx  +  cos  9mc z)  •  {pi  cos  </>x  +  p{  sin  <j>y  +  ztz) 
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=  -kpi  sin  9inc  cos  0  -  kzi  cos  9mc  (55) 

where  as  in  Section  2.3  we  have  let  pi  —  p(t{ )  and  z%  —  z(ti).  Hence,  recalling  that  it  has 
been  assumed  that  0inc  =  0,  we  then  obtain  from  (53) 

2t r 

yte  =  keJ^iCos0^coseincsinVij  d(p  cos  (fjep&Pi sin  °inc  cos 

0 


i/I  _ 

-  sin  9inc  cos  Vi  I  d<f>ej(kpiSin9inCcos4,-n,l’)  J 
0 

2n 

V%  =  -  cos 9inckejkziCOseinc  J d<psin(f>eKkpi8inei 

0 
27T 

Vjf  =  kePkZiCOs6incsmVij  d<P  sin  4>ej{kp'sindinCcos4,-n,p) 


$inc  cos  <j>—n<f>) 


and 


27T 

Vff  =  h(jkz'cos0inc  J  d(f) cos0ej(fcpiSi 


sin  0tnc  cos  <j>~n<j>) 


(56a) 

(56b) 

(56c) 

(56d) 


where  we  have  changed  the  limits  of  the  (f)  integration  to  0  to  2tt.  Using  the  integral  repre¬ 
sentation  for  the  Bessel  function 


_ n  2n 

JJx)  =  3—j 


(57) 


it  is  simple  to  obtain  the  integral  representations 


and 


so  that 


and 


2tt 

J  sin  =  -7rf[Jn+i(:r)  +  Jn- i(x)} 

0 

27T 

J  coscjxj^f-^dt  =  7T jn+1[Jn+l(x)  -  Jn-l(x)] 

0 

V"  =  Ttjn+,k  [cos9i“sinui(Jn+1  -  J„-i)  +  2jsin9''“cost',J„]  e?tn 

Vff  =  «»«*“( 

V„"f  =  -»/**» in^(/„+.  +  Jn-i)^^’ 

V*?  =  -  Jn- 


cos  0*nc 


(58a) 


(58b) 


(59a) 

(59b) 

(59c) 

(59d) 
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where  we  have  let 

Jn  =  Jn(kpi  sin#”*). 

(59e) 

Using  the  Bessel  function  relation 

J-n(z)  =  (-!)"./„(*) 

(60) 

it  is  easy  to  show  from  (59)  that 

£ 

H 

**• 

(61a) 

v —ni  vni 

(61b) 

yty  _  _ 

v  —ni  v  ni 

(61c) 

and 

\ r<t><P  t r<M> 

v  —ni  vni  * 

(61d) 

2.6  Calculation  of  the  Current 

To  calculate  the  current  on  the  BOR  surface  we  refer  to  (19)  and  (24)  and  write 

K «(*,*)=  £  t**  £  [*&*(*){((,*)  +  *&(*)*(*)]  (62) 

n——N  i 


where  q  =  6  ( TM )  or  <j>  (TE)  indicates  the  polarization  of  the  incident  electric  field  given 
by  (16)  or  (17),  and  (K%K%[)  are  the  elements  of  the  vectors  (Ktq  ,K£q) ,  obtained  as  the 
solution  of  the  matrix  equation  (see  (25)) 


\K\ 

IS?]  ' 

'K?  ' 

'  \r*<i 

v  n 

\zf\ 

[ZH]  J 

kh 

\r<k 

n 

n  —  0,  ±1,  ±2,  •  •  • ,  iN. 


(63) 


Letting  p  =  \pi(t),P2(t),  •  •  -PN(t)]  be  the  row  vector  of  the  pulse  basis  functions  pi(t)  defined 
by  (23)  (we  use  the  tilde  to  indicate  the  transpose  of  a  column  vector),  (62)  can  be  written 
as 


K’M)=  £; 


(64) 


n=-N 


From  (63),  the  relations  (57)  for  the  Z_n  matrices  in  terms  of  the  Zn  matrices,  and  the 
relations  (61)  for  the  V-n  vectors  in  terms  of  the  Vn  vectors,  it  is  simple  to  derive  the 
relations 


and 


_ J£t0 

(65a) 

1 

II 

(65b) 

=  -K* 

(65c) 

B 

3^ 

II 

B 

(65d) 
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(Note  that  (65b)  and  (65c)  imply  that  Kf  =  0  and  K?  =  0.)  Substituting  (65)  in  (64)  then 
yields  the  expressions  for  the  currents  induced  on  the  surface  of  the  BOR  by  a  TM  and  TE 
linearly  polarized  incident  plane  wave,  respectively: 

Ke(t , 0)  =  pK$t(t, <t>)  +  2  Y,  [(PKn) cos n&  +  jipK-n) sin n(t>H<f>)\  (66a) 

n= 1 

and  N 

K +(t,  <f>)  =  pKf4>  +  2  £)  [(pKj?) sin  0)  +  (^n0) cos  ™M>(<£)]  •  (66b) 

n=l 

2.7  Calculation  of  the  Far  Scattered  Field 

In  the  fax-field  region  the  9  and  <fi  components  of  the  scattered  electric  field  can  be  obtained 
from  the  vector  potential  A  [10,  p.  281] 

Esc(r)  r~°°  -ju;A(r).  (67) 


Then  with  (3)  and  (5), 

jfc|  r— r'| 

Hi)  -3kz 0 /[K’(r')  •  p]4^  Ir-r,  |dg'  (68) 

s  1 

*  xfar  lfar  • 

where  p  :  9  or  0,  indicates  the  far  scattered  field  component;  p  :  0  or  <p  ,  is  tne 
corresponding  unit  vector  (see  Fig.  2);  q  :  9  or  0,  indicates  the  polarization,  9  (TM)  or 
<f>  (TE)  of  the  incident  electric  field;  K«  is  given  by  (64),  r  =  (r,  0/ar,  (j>far)  is  the  far-field 
observation  point,  and  r'  =  (p' ,  <//,  z')  is  the  source  point  on  the  BOR  surface.  In  the  far 
field,  |  r  -  r'  |  in  the  phase  of  the  Green’s  function  can  be  approximated  by 

Ir-r'  f  ~°°  r  —  p'  sin  9far  cos(0'  -  4>far)  -  z'  cos  6far  (69) 

so  that 

Esc™{ r)  r  ~°°  —jkZo^——  J [K9(r')  •  p\eS*ff^^^-^ar)etkl'COB$,ardS'.  (70) 

s 

The  relations  (54)  can  be  used  directly  to  obtain 

t*  •  0^ar  —  cos  0far  sin  v(t')  cos (<f/  —  <p*ar)  —  sin  9^ar  cos  v(t')  (71a) 

<j>  ■  $^ar  =  —  cos  9^ar  sin(</>/  —  <f>far)  (71b) 

t#  •  (f/ar  =  sin  v{t')  sin  (ft  —  <j>iar)  (71c) 

and  , 

<j>  •  <j>  ar  =  cos(4>'  —  4>far)-  (71d) 
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Letting  the  elements  of  the  vectors  Rjf  and  be  defined  by 


J pi(t')( t#  •  p)^*^ sin efar cos(<t,,-<t,far)+z' c™ °lar] -<t>far) dS' 
s 

and 

R%  =  kJ  Pi(t')(4 >'  •  p)e>k[f>' sin  6far  <**(*' -<t>/ar)+*'  “s  efar] 

s 

with  p  —  $  or  (f>,  (70)  together  with  (64)  can  be  written  as 

Esc^{r)  r  ~°°  ~jkZ°e~3kr  £  eJ'^/or(R^K^  + 

Airkr  „=_# 

where,  using  (72)  and  the  2n  periodicity  of  the  4>'  integration 


(72a) 

(72b) 


(73) 


2-k 

R%  =kj  dt'p'piit')^2'™5  efaT  [  cos  efar  sin  v(t')  J  d<j>'  cos  4>e?{kp' sin  0far  003  *'+n*') 

o 

277 

-sin«/arcost>(f')  J  ^V(V 


(74a) 


27T 

R%  =  -  cos  efark  J  dt'p'piit ')e?kz,c°sOfaT  J  d<j>'  sin  ^V'(^'sin0/arc^W)  (74b) 


2tv 

R%  =  kJ  dt'p'pi  (t')e?kz>  cos  efar  smv(t')  J  d<f>’ sinfiepW™6'"008*'^  (74c) 

o 

and 

2n 

Rtt  =  kj  dt,p'pi{t')ejkz'cose,ar  j  d(j)'  cos  4>'(j{kp' sin  efar  005  ^'+n*') .  (74d) 

J  o 

Comparing  (74)  with  (56),  replacing  6inc  of  (56)  by  0far  and  —n  of  (56)  with  n,  and  using 
the  definition  of  the  pulse  basis  functions,  we  then  obtain 


R%  =  7 Tjn+1kdiPi  [cos  efar  sin  7 H(Jn+i  -  <4-1 )  +  2 j  sin  9far  cos  ViJn\  eJ'fcz‘cose/°r  (75a) 

R%  =  —njnkdiPi  cos  0^ar(  Jn+1  +  Jn-x)^^008^,  (75b) 

=  njnkdipisinvi(Jn+1  +  Jn-1)eJ'*‘<C0S*/or  (75c) 

and 

fljf  =  njn+1kdipi(Jn+i  -  J„-i)e’l““x’e'"  (75d) 

where 

<4  =  <4(/cPi  sin  0'ar).  (75e) 
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The  relations  (61)  apply  to  the  Rp f  as  well,  so  that 

rytO  _  -ntB 


r><t>Q  _  _ 

■**'— m  iLm 

/?!,•  =  -i$ 


and 


iy4>tfr  _  c 

ni  lhni  * 

Then,  substituting  (76)  and  (65)  in  (73)  it  is  straightforward  to  obtain 


Esc’"(r) 


r— *  OO 

ryj 


jkZoe 


-jkr 


47T  kr 


+  2  £  (R?K?  +  R*Kj*)  cos(n^,“r) 

n—1 


E 


sc,<P6(r)  r^oo 


kZ0e~ik\  " 


£sc'^(r) 


1 - KX) 


,  ,  2£(it?Kl?+R*Kf)sm(n/‘”') 

47T/cr  n==1  v 

1^2  £  (R?K?  +  RfKj*)  sta(n/-) 


n=l 


and 


gsc,<l>4>^ 


r— ►oo 


jkZ0e 


-jkr 


+  2  f;  (r*K*  +  R$*KJ*)  cos(n0/ar) 

n=l 


47t  kr 

The  radar  cross  section  cr  is  defined  as 

apq  —  lim  47rr 


2 1  Esc'pq  |2 

ginc,g  2 


(76a) 

(76b) 

(76c) 

(76d) 

(77a) 

(77b) 

(77c) 

(77d) 


(78) 


where  p  :  6  or  <f>,  denotes  the  component  of  the  far  scattered  field,  and  q  :  9  (TM)  or  cj)  (TE) , 
indicates  the  polarization  of  the  incident  electric  field. 


2.8  Choice  of  the  Number  of  Fourier  Modes  for  Expansion  and 
Testing  Functions. 

The  4>  dependence  of  the  current  given  by  (19)  is  expressed  as  a  summation  from  —N  to  N 
of  the  Fourier  modes  ejn<t>.  The  value  of  N  can  be  set  equal  to  the  number  of  Fourier  modes 
sufficient  to  represent,  to  the  desired  accuracy,  the  <j>  variation  of  the  tangential  component 
of  the  incident  electric  field  on  the  surface  of  the  BOR.  Let  a  be  the  largest  value  of  p  of  a 
point  (p,  z)  on  the  generating  curve  of  the  BOR.  Then  from  (15),  (16),  and  (17)  it  can  be 
seen  that  the  <f>  variation  of  the  incident  field  along  the  circle  on  the  BOR  corresponding  to 
the  point  (a,  z)  is  given  by 

=  cos<^'fcasin®<nCcos*  or  =  sin<^'fcosin*<nCcos*.  (79) 
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For  axial  incidence  6inc  =  0  and  f{4>)  equals  cos  (f)  or  sin  4>  so  that  only  the  e±J^  modes  are 
needed.  For  oblique  incidence,  we  can  express  f(<p)  as  the  Fourier  series 

m-  Ew-  (80a) 

-N 


with 

*  2?T 

ft.  =  2 -  / (80b) 
0 

from  which,  with  (79)  and  (58), 

\cn\=\\  Jn+i(kasm0inc)  ±  Jn^(ka  sin  0inc)  |  .  (81) 

As  n  increases  beyond  fca  sin  dtnc,  |  Jn(ka  sin  0inc)  |  decreases  rapidly  and  it  suffices  to  choose 

N  =  I  +  M  (82) 


where 

I  —  Int[(l  +  a)fcasin0mc],  0<aCl 
and  M  is  the  smallest  integer  for  which 

J^(ka  sin  9mc) 

Ji(ka  sin  6inc)  ~  6 


(83) 

(84) 


with  e  a  small  positive  number  depending  on  the  desired  accuracy.  If  the  value 
by  (82)  is  plotted  as  a  function  of  Int[/ca  sin  #wc]  the  plot  is  found  to  be  almost 
e  =  0.005,  for  example, 


N  «  Int[l.O4fcasin0inc]  +  7. 


of  N  given 
linear.  For 

(85) 


A  similar  expression 


N  ~  Int[/c+(asin  6)  +  A]  (86) 

with  k+  denoting  a  value  a  few  percent  larger  than  k  was  obtained  by  Yaghjian  [11]  for  the 
reciprocal  problem  of  estimating  the  number  of  angular  modes  needed  to  represent  the  far 
field  of  a  radiator  in  the  6  direction. 


2.9  Approximation  of  Derivatives 

To  approximate  the  derivatives  in  (36),  (43),  and  (48)  we  employ  a  quadratic  Lagrangian 
interpolation  polynomial  [12].  Recall  that  in  (24)  Kln{t)  and  K%(t)  were  approximated  using 
pulse  basis  functions.  The  values  of  pK^{t)  and  (1  /p)K%(t)  at  the  discrete  values  t  =  tj,j  = 
1,  •  •  • ,  M  must  then  serve  as  the  basis  for  approximating  the  derivatives.  To  do  this  we  fit 
quadratic  polynomials  g*(t)  and  q^it)  through  three  consecutive  values  of  pKl^  and  (1  /p)K% 
respectively: 


<f(t)  =  Pi-iKfo-i) 


{t  tj)(t  tj+i) 

{tj~~  i _  ^7X^7— 1  ”  1) 


(t  tj+i) 

{pj  ~~  “  tj+ 1) 
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These  approximations  are  of  the  form  (37)  and  (38). 


(tj-i  —  tj)(tj~i  —  tj+i)  ( tj+i  tj-i)(tj+i  tj) 

^Pj+i^n(tj+i)  (88c) 


3  NUMERICAL  RESULTS 

The  analysis  presented  in  Section  2  was  implemented  in  a  FORTRAN  computer  program 
which  was  then  used  to  obtain  numerical  results  for  several  different  BOR’s.  In  this  section 
we  show  some  representative  results  using  the  LEFIE  to  calculate  scattering  from  a  prolate 
spheroid.  The  geometry  of  the  spheroid  is  shown  in  Figure  4.  The  semi-major  axis  of  the 
spheroid  is  a  and  the  semi-minor  axis  is  b.  For  the  calculations  we  show,  ka  =  20  and 
kb  =  10.  The  TM  illuminating  plane  wave  makes  an  angle  of  45°  with  the  major  axis  of  the 
spheroid.  For  these  values  of  ka  and  kb  no  spurious  resonances  are  encountered.  (See  end  of 
Section  2.2.) 

The  solid  curve  in  Figure  5  shows  the  E-plane  pattern  calculated  with  the  combined 
field  integral  equation  (CFIE)  implemented  in  the  computer  code  CICERO  [13]  using  a 
discretization  of  the  spheroid  generating  curve  of  40  points/ A.  The  dotted  curve  in  Figure 
5  is  the  pattern  obtained  by  solving  the  conventional  EFIE  with  the  Galerkin  form  of  the 
method  of  moments  and  overlapping  triangle  basis  functions  with  a  point  density  of  20 
points/ A.  The  dot-dashed  curve  and  the  dashed  curve  in  Figure  5  are  the  patterns  obtained 
solving 
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Figure  5:  E-plane  pattern  of  a  prolate  spheroid  with  ka  —  20  and  kb  —  10  illuminated  by  a 
TM  plane  wave  incident  at  an  angle  of  45°  with  the  major  axis,  as  calculated  by  the  CFIE 
at  a  density  of  40  points/A,  by  the  EFIE  solved  with  overlapping  triangle  basis  functions 
at  a  density  of  20  points/A,  and  by  the  EFIE  solved  with  pulse  basis  functions  and  point 
matching  with  a  density  of  20  points/ A  and  80  points/ A. 
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the  conventional  EFIE  with  pulse  basis  functions  and  point  matching  at  point  densities  of 
20  points/A  and  80  points/A,  respectively.  We  note  that  the  CFIE  and  EFIE  solved  with 
overlapping  triangle  basis  functions  yield  results  that  agree  closely,  but  the  patterns  obtained 
with  pulse  basis  functions  and  point  matching  are  highly  inaccurate.  Moreover,  these  pulse 
basis  function  patterns  are  not  significantly  improved  upon  by  considerably  increasing  the 
point  density. 

In  contrast  with  the  very  poor  results  obtained  when  the  conventional  EFIE  is  solved 
using  pulse  basis  functions  and  point  matching,  in  Figure  6  we  show  the  spheroid  patterns 
obtained  solving  the  conventional  EFIE  with  the  method  of  moments  and  overlapping  tri¬ 
angle  basis  functions  with  a  density  of  20  points/A,  and  the  patterns  obtained  by  solving 
the  LEFIE  with  pulse  basis  functions  and  point  matching  at  densities  of  40  pointsA  and  20 
points/A.  The  LEFIE  solved  at  a  density  of  40  points/A  yields  a  pattern  quite  close  to  the 
conventional  EFIE  pattern  obtained  at  a  density  of  20  points/ A.  The  LEFIE  solved  with  a 
density  of  20  points/A,  however,  has  some  pattern  errors  of  approximately  1.5  dB.  The  point 
density  required  to  obtain  results  with  the  LEFIE  comparable  in  accuracy  to  those  obtained 
with  the  conventional  EFIE  can  be  improved  upon  considerably  by  the  stratagem  of  using 
a  higher  point  density  for  only  a  very  small  region  —  say  A/4  -  in  the  vicinity  of  the  ends  of 
the  spheroid  generating  curve,  and  a  low  point  density  elsewhere  as  shown  in  Figure  7.  In 
the  LEFIE  pattern  shown  in  Figure  7  a  density  of  80  points/A  was  used  close  to  the  ends  of 
the  spheroid  generating  curve  (i.e.,  the  poles  of  the  spheroid)  and  a  density  of  20  points/A 
was  used  elsewhere.  The  reason  why  a  high  point  density  may  be  required  at  the  ends  of  the 
BOR  generating  curve  is  that  the  finite  difference  approximation  of  the  derivatives  used  in 
solving  the  LEFIE  is  less  accurate  at  the  beginning  and  end  of  the  generating  curve  because 
one-sided  derivative  approximations  must  be  used  there.  Using  a  higher  point  density  at  the 
ends  of  the  generating  curve  compensates  for  the  use  of  one-sided  finite  difference  derivative 
approximations. 

4  SUMMARY 

Unlike  the  magnetic-field  integral  equation,  the  conventional  form  of  the  electric-field  integral 
equation  cannot  be  solved  accurately  using  pulse-basis  functions  and  point  matching.  It 
can  be  demonstrated  that  it  is  the  highly  singular  kernel  of  the  EFIE,  rather  than  the 
derivatives  of  the  current,  that  precludes  the  use  of  the  pulse-basis-function,  point-matching 
MOM.  A  new  form  of  the  EFIE  has  been  derived  whose  kernel  has  no  greater  singularity 
than  the  free-space  Green’s  function.  This  new  low-order  singularity  form  of  the  EFIE,  the 
LEFIE,  has  been  solved  for  a  perfectly  electrically  conducting  body  of  revolution  using  the 
pulse-basis  function  point-matching  MOM,  and  a  computer  program  has  been  written  to 
implement  the  solution.  Derivatives  of  the  current  axe  approximated  with  finite  differences 
using  a  quadratic  Lagrangian  interpolation  polynomial.  This  simple  solution  of  the  LEFIE 
is  contingent,  however,  on  the  vanishing  of  a  linear  integration  term  that  appears  when  the 
original  EFIE  is  transformed  to  obtain  the  LEFIE.  The  vanishing  of  this  linear  integration 
term  restricts  the  applicability  of  the  LEFIE  to  smooth  closed  scatterers  without  tips,  edges, 
or  any  other  features  that  might  result  in  a  discontinuous  or  unbounded  surface  charge 
density,  its  surface  derivative,  and  the  sum  of  the  reciprocals  of  the  principal  radii  of 
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Figure  6:  El-plane  pattern  of  a  prolate  spheroid  with  ka  —  20  and  kb  —  10  illuminated 
by  a  TM  plane  wave  incident  at  an  angle  of  45°  with  the  major  axis,  as  calculated  by  the 
EFIE  solved  with  overlapping  triangle  basis  functions  at  a  density  of  20  points/ A,  and  by 
the  LEFIE  solved  with  pulse  basis  functions  and  point  matching  at  a  density  of  40  points/A 
and  20  points/ A. 
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Figure  7:  E-plane  pattern  of  a  prolate  spheroid  with  ka  =  20  and  kb  —  10  illuminated 
by  a  TM  plane  wave  incident  at  an  angle  of  45°  with  the  major  axis,  as  calculated  by  the 
EFIE  solved  with  overlapping  triangle  basis  functions  at  a  density  of  20  points/A,  by  the 
LEFIE  solved  with  pulse  basis  functions  and  point  matching  at  a  density  of  80  points/A  in 
an  interval  of  A/4  at  either  end  of  the  spheroid  generating  curve  and  20  points/A  elsewhere, 
and  by  the  LEFIE  at  a  uniform  density  of  20  points/A. 
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curvature,  or  a  discontinuous  unit  normal  vector  to  the  surface.  Bistatic  RCS  calculations 
performed  for  a  prolate  spheroid  demonstrate  that  results  comparable  in  accuracy  to  the 
conventional  EFIE  can  be  obtained  with  the  LEFIE  using  pulse-basis  functions  and  point 
matching  provided  that  a  higher  density  of  points  is  used  close  to  the  ends  of  the  BOR 
generating  curve  to  compensate  for  the  use  of  one-sided  finite  difference  approximations  of 
the  first  and  second  derivatives  of  the  current. 


References 

[1]  A.W.  Maue,  “On  the  formulation  of  a  general  scattering  problem  by  means  of  an  integral 
equation,”  Z.  Phys.,  vol.  126,  pp.  601-618,  1949. 

[2]  J.  Van  Bladel,  Electromagnetic  Fields,  McGraw-Hill,  New  York,  1964. 

[3]  A.D.  Yaghjian,  “Augmented  electric-  and  magnetic-field  integral  equations,”  Radio  Sci¬ 
ence,  vol.  16,  987-1001,  November-December  1981. 

[4]  J.R.  Mautz  and  R.F.  Harrington,  “H-field,  E-field,  and  combined-field  solutions  for 
conducting  bodies  of  revolution’”  Arch.  Math.  Ubertragungtech.,  vol.  32,  pp.  157-164, 
1978. 

[5]  A.D.  Yaghjian  and  M.B.  Woodworth,  “Derivation,  application  and  conjugate  gradient 
solution  of  the  dual-surface  integral  equations  for  three-dimensional,  multiwavelength 
perfect  conductors,”  in  PIER  5,  Progress  in  Electromagnetics  Research:  Applications 
of  the  Conjugate  Gradient  Method  to  Electromagnetics  and  Signal  Analysis,  T.  Sarkar, 
ed.,  New  York:  Elsevier,  1991,  103-130. 

[6]  M.B.  Woodworth  and  A.D.  Yaghjian,  “Multiwavelength  three-dimensional  scattering 
with  dual-surface  integral  equations,”  J.  Opt.  Soc.  Am.  A,  vol.  11,  pp.  1399-1413,  April 
1994. 

[7]  A.F  Peterson,  S.L.  Ray,  and  R.  Mittra,  Computational  Methods  for  Electromagnetics, 
New  York:  IEEE  Press,  1998. 

[8]  R.A.  Shore  and  A.D.  Yaghjian,  “ Dual  Surface  Electric  Field  Equation”,  Air  Force  Re¬ 
search  Laboratory  Report  No.  AFRL-SN-HS-TR- 2001-013,  2001. 

[9]  R.A.  Shore  and  A.D.  Yaghjian,  “Dual-Surface  Integral  Equations  in  Electromagnetic 
Scattering,”  Proceedings  International  Union  of  Radio  Science  (URSI)  27th  General 
Assembly,  Maastricht,  The  Netherlands,  August  2002. 

[10]  C.A.  Balanis,  Advanced  Engineering  Electromagnetics,  New  York:  John  Wiley,  1989. 

[11]  A.D.  Yaghjian,  Near  Field  Antenna  Measurements  on  a  Cylindrical  Surface:  A  Source 
Scattering  Matrix  Formulation,  NBS  Technical  Note  696,  1977. 

[12]  C.E.  Pearson,  Numerical  Methods  in  Engineering  and  Science,  New  York:  Van  Nostrand 
Reinhold,  1986. 

[13]  J.M.  Putnam  and  L.N.  Medgyesi-Mitschang,  Combined  Field  Integral  Equation  Formu¬ 
lation  for  Axially  Inhomogeneous  Bodies  of  Revolution,  McDonnell  Douglas  Research 
Laboratories  MDC  Report  No.  QA003,  December  1987. 


31 


